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Determinaiion of the elemental composition of cosmic dust particles by means of an 
impact ionization time-of -flight mass spectrometer has been investigated (d several in- 
stitutions. In most configurations, the mstrument supplies the identity of ion groups of 
both target and particle materials extracted from the impact plasma and the number of 
ions contained in each group. Experiments have shown that the frewtional ionization 
of a given species is not constant with impact velocity nor is thefra^ional ionization the 
same for different kinds of atoms. A model of the impact ionization effect developed at 
TRW involves an equilibrium plasma condition with the consequence that the fractional 
ionization for an arbitrary atomic species can be specified by the Saha equation if the 
plasma volume (V) and temperature (T) are known. It follows that T can be determined 
by taking the ratio of the Saha equations for two elements present in the target in known 
concentration. (Taking the ratio negates the requirement of knowing V.) Given T, the 
procedure can be reversed to yield the relative abundance of elements contained in the 
impacting particle. To test the model, a PbZrOsrPbTiOi target wasbombarded with high 
velocity Fe, MoB, and NiAl particles and the number of Pb, Ti, and Zr ions was de- 
termined in a time-of-flight mass spectrometer. For each event, the relative abundance 
of Ti to Pb teas taken as known (from electron microprobe analysis) and T was deter- 
mined from the Ti-Pb measurement. The Zr to Pb ratio was found to be in good agree- 
ment with the microprobe analysis (0.38 calculated mean value compared to 0.34 actual). 
The result ivas valid for all particle materials and for a velocity range 17 <v<47 hm/s. 
T ranged from 3300 to 11 500° K and was only mildly velocity dependent. 


T he cosmic dust analyzer under development 
at TRW may be described as an impact 
ionization time-of-fligbt mass spectrometer. As 
implied by this description, the instrument 
utilizes ion time-of-flight techniques to separate 
the microplasma formed by the impact of a high 
velocity particle upon a solid surface into its 
constituent parts. The plasma contains ions of 
both particle and target material, which may be 
extracted from the impact plasma by the applica- 


tion of an electric field to the impact surface. Tliis 
is most easily accomplished by placing a potential 
difference between the impact surface and a high 
transmittance grid located in close proximity to 
the surface and parallel to it. All singly ionized 
atoms reach the plane of the grid with the same 
kinetic energy (given by the product of the unit 
electronic charge and the applied voltage) ; how- 
ever, the velocity of a particular ion at this point 
depends upon its mass. Thus, the atomic mass of 
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the ions may be determined by measuring their 
transit time from the grid to an ion collector 
located some distance away. 

First, the formation of the plasma and the 
extraction of ions from it may be considered 
instantaneous, with the conscQuence that all ions 
of a given mass reach the collector at the same 
time. Actually, it has been observed that the ions 
may possess an appreciable thermal (random) 
velocity which adds to the electrically imparted 
velocity. This causes a small spread in arrival 
times of ions of a given chargo-to-mass ratio. 
Nevertheless, by proper selection of accelerating 
voltage and drift distance, the ions can be segre- 
gated into distinct groups according to their mass. 
With the accelerating voltage and flight distance 
specified, measurement of the transit time is 
sufficient to determine the mass of each ion group. 
Determination of the atomic mass identifies the 
corresponding element. 

The foregoing defines the basic instrumental 
requirements of the cosmic dust analyzer. The key 
measurements to be made are the transit time of 
the various ion groups from the impact surface to 
the ion collector and the total charge contained 
within each group. The instant of impact is 
signaled by the emission of positive charge from 
the impact (target) surface, resulting in a negative 
signal. This signal may be used to initiate the ion 
time-of-flight measurement. A convenient means 
of measuring the total charge within each ion 
group is to integrate the ion current at the collector. 
In this case the ion collector signal has a stairstep 
appearance; the position in time of a step corre- 
sponds to the time of arrival of a given ion group 
and its amplitude is proportional to the charge 
contained within the ion group. The number of 
ions in each group is obtained by dividing the 
total charge by the unit electron charge. 

According to its description, the instrument 
supplies direct measurements of the number of 
ions in each group reaching the collector and 
provides the timing information necessary to 
identify the species of the ions. These data alone 
would be extremely useful in evaluating the gross 
characteristics of cosmic dust particles. However, 
a more sophisticated obj (active is that of quantita- 
tively determining the relative abundances of the 
neutral elements contained in the particles. As 
might be expected intuitively and as has been 


verified experimentally, the fractional degree of 
ionization of all atomic species is not the same. 
Moreover, the relative ionizing efficiencies change 
with impact velocity. Hence, additional informa- 
tion is needed in order to process the raw data to 
obtain the desired final result. A method has been 
conceived to obtain that information; this paper 
will report the results of laboratory experiments 
which indicate that the method is a valid approach 
to the problem. 

At the heart of the method is a theoretical model 
of the impact ionization process that was origi- 
nally developed because it successfully explained 
various phenomena observed in earlier impact 
experiments. Briefly, the model assumes that the 
impact produces a plasma in thermal equilibrium. 
If this is the case and if the plasma temperature 
and volume can be found, the degree of ionization 
of any atomic species in the plasma can be com- 
puted provided that the species has a known 
ionization potential. The requirement to specify 
the plasma volume can be shown to disappear if 
we arc content to determine relative abundances, 
leaving only the plasma temperature to be meas- 
ured, If the impact target contains two elements in 
known concentration, the model provides equa- 
tions by which measurement of the relative ion 
signals of these two elements can be reduced to 
plasma temperature. Given the temperature, the 
model supplies further theoretical relationships for 
processing the data to yield the relative atomic 
abundances of the elemental constituents in the 
impacting particle. 

This model of impact ionization and the 
theoretical relationships derived therefrom should 
be experimentally verified. The best approach to 
this task appears to be to employ test particles 
and an impact target that are each composed of 
two elements in known relative concentrations, 
utilizing the theory to be described in the section 
of this paper on the impact ionization model to 
independently determine the plasma temperature 
from each of the two-element combinations. If 
the two values of plasma temperature so deter- 
mined should agree for each hypervclocity impact 
event, then the model may be considered to have 
been verified. Unfortunately, we have been unable 
to locate a source of particles with sufficiently 
homogeneous properties to successfully complete 
this type of experiment. 
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We have therefore developed an alternative 
approach that seems to be almost as satisfactory. 
The alternate involves the use of a target material 
containing three major elements Mith known 
relative abundances. Tw’o of these arc treated as 
control elements, and the relative ion abundances 
of the control species are observed and used to 
compute the plasma temperature as previously 
described. This plasma temperature and the 
magnitude of the ion signal from the tliird clement 
are then combined to compute the relative atomic 
abundance of the latter. The computed value can 
of course be compared directly to the known 
abundance of the third material; good agreement 
indicates that the model is essentially correct. 
In comparison with the postulated ^^best ap- 
proach, the alternate suffers by not directly 
involving the properties of the impacting particle 
in the calculations, and the verification provided 
thereby assumes that the target and particle 
materials become intimately mixed in the impact 
plasma. This assumption requires its own experi- 
mental test. Even in the absence of such a test, 
hownyer, any agreement found from the alternate 
experiment is strong evidence that the model is a 
reasonably good representation of the impact 
ionization phenomenon. 

We will describe herein an experiment in which 
the alternative approach wns employed and in 
which data w^'cre acquired for a number of impact 
events involving different particle compositions. 
As we hope to show^, these data appear to possess 
a sufficiently good agreement with the known 
target composition that the desired verification 
can be said to have been achieved. 

IMPACT IONIZATION MODEL 

It w^as first shown experimentally by Friich- 
tenicht and Slattery (1963) that free electric 
charge, both positive and negative, is produced at 
the point of impact of a high velocity micro- 
particle on a solid surface. More recent experi- 
ments reported by Hansen (1968) and by Auer 
and Sitte (1968) have shown that the positive 
charge consists of ions of both the impacting 
particle and target materials and that the negative 
charge carriers are electrons. Unpublished data 
obtained at TRW have demonstrated this effect 
in considerable detail for a large number of 


particle-target material combinations. For low^ 
impact velocities, both Hansen and Auer and 
Sitte found that easily ionized impurities (sodium 
and potassium) accounted for the majority of the 
emitted positive charge even though the impacting 
particles wrere ^^98 percent iron. At higher impact 
velocities, how^ever, iron became the dominant 
ionic species, and ions of other high ionization 
potential elements made an appearance. Particu- 
larly significant w^as the detection of ions from 
refractorj^ metal targets. 

Based on these observations, Hansen (1968) 
advanced a theoretical model to describe the 
impact ionization effect. In this model, it is 
assumed that a small, high density plasma is 
formed in the immediate vicinity of the impact 
site. The plasma expands due to its internal 
pressure; however, the collision frequency of 
atoms and electrons within the plasma is very 
high (on the order of 10^^ per second initially), 
with the consequence that the plasma tempera- 
ture is spatially uniform (i.e., that the plasma 
temperature is nearly the same throughout the 
plasma volume) but decreases a^ the plasma 
expands. At some point the plasma becomes 
collisionless as the density decreases. To the 
extent that the plasma may be considered to be in 
thermal equilibrium immediately prior to be- 
coming collisionless, the Hansen model described 
in the following paragraphs should apply. 

When a plasma in thermal equilibrium at 
absolute temperature T contains several species of 
atoms and ions, the degree of ionization of each 
species s whose ionization energy is Ej^ is related 
to T according to Saha^s equation (cf. Sutton and 
Sherman, 1965) : 

= expi-Ei‘/kT) (1) 

where Uej and Un are the numbers per unit 
volume of electrons, ions of species s, and neutral 
atoms of species s, respectively; C is a constant, 
Ui,n’‘ is the ratio in which Ui* and Wn* are 

respectively the internal partition functions of 
ions and neutrals of species s; and h is Boltzmann s 
constant. The Saha equation show's that the 
fractional ionization of a given species of atom 
depends on the ratio of the ionization energy of 
that species to the mean thermal energy of the 
plasma. For reasonably low^ plasma temperatures, 
only a small fraction of the atoms are ionized 
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(except for those with very low ionization poten- 
tials) . In order to use the impact ionization effect 
to reliably determine the relative abundance of 
the elements in cosmic dust and meteoroids, a 
measure of the plasma temperature must be 
obtained. 

It would at first appear that the number density 
of ions, electrons, and atoms must be specified. 
However, from the argument given above, the 
plasma rapidly expands to a collisionless state. 
Thus, the total number of ions is constant at later 
times during the expansion process. Since all of 
the ions regardless of species are extracted from 
the same volume, a direct measurement of Ue, 
7?n*, and is not required to determine the 
relative abundance of the different kinds of atoms. 
If iVc, Nn% and iV,* are, respectively, the total 
numbers of electrons, atoms, and ions in the 
plasma, and if V is the volume of the plasma at any 
time prior to its expansion to a collisionless state, 
then Ue = Ne/V, = and ni* = Ni^/V, 

Clearly, Ueni^lun* can be replaced in equation (1) 
by ’NJSixj^nY, As we will shortly sec, F can be 
ignored in the typa of experiment that we perform. 
Our experimental techniques provide a direct 
measurement of for each 5 and of Ne. (If all 
are measured, can be obtained as 
since the total number of electrons is equal to that 
of ions; however, an independent measurement of 
A e is readily obtained and serves as a useful check 
of the data.) The quantity of principal interest to 
be extracted from the experiments is 

Consider first an idealized case in which the 
plasma contains two species from the target and 
one from the particle; let these species be denoted 
by superscripts ( 1) , (2) , and (3) , and assume that 
the identities of (1) and (2), their ionization 
energies and and their partition function 
ratios and are all kno\\m. If equation 
(1) is written for species (1) and for species (2) 
and if the former result is divided by the latter, 
one finds that 

(A^<^>/A^c 2>) (A^„(2)/^„a>) 

=/'(T) exp £(Ej^^^-~Ej(^))/kT:\ (2) 

in which /'(T) (The partition 

functions are functions of temperature; the func- 
tional dependence is assumed to be known.) It is 
almost always true that Ar„*»A^*, and then the 
ratio is approximately equal to the 


relative concentration of the two kinds of atoms 
in the target material. Let A denote this ratio, 
whose value is known. Now, suppose that species 
(1) and (2) have been selected such that 
EY^K Writing = sub- 

stituting into equation (2), we obtain 

exp( — A£'i<^'2y^iT) (3) 

from which 

kT=AEi^^>^^]n C/'(T)A(AVWi<'>)]}-^ (4) 

According to equation (4), the plasma tempera- 
ture can be found from a measurement of the 
number of ions of the two control elements in the 
target. (Remember that/'(T) is a kno'^ra func- 
tion.) If, as often happens, f'(T) is essentially 
constant over a range of temperatures, the com- 
putation of T from the measurement is 

quite simple; if f'(T) is not constant, iterative 
solutions can be employed. 

It follows at once that the procedure employed 
to derive equation (2) can be applied for species 
(1) and (3), with the result: 

Xexp[(iJ,(3)-£J,a))/A:r] ( 5 ) 

w^here g'{T) - Thus, if Wi^VA^/D has 

also been measured, the relative abundance 
can be computed. The procedure is 
readily generalized to a situation in w^hich the 
particle contains a multiplicity of elements and the 
relative abundance of each is computed in turn 
from measurements of the corresponding ion 
ratios. This is the way in which it is intended to 
reduce data from the actual cosmic dust analyzer 
experiment. 

We saw that one way of checking the validity 
of these equations in this application w’ould be to 
employ not only a target containing two control 
species (1) and (2), but also a test particle con- 
taining two other species %vhose relative concen- 
trations in the particle are known; these are 
species (3) and (4) . If the impact ionization model 
is valid, an equation similar to equation (4) can 
be WTitten for the latter species, replacing 
A^i< 2 )/Arxi) by W,wyv,-(3), f(T) by A'(T) = 
and AEY^’^^ by AEi^^’^K Then a 
measurement of would give an in- 

dependent determination of T that could be 
compared with that found from A^^ 2 )/^YXd. agree- 
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ment would strongly imply that the predictions of 
the model are correct. 

The alternative check, used as the basis of our 
experiment for the reasons stated earlier, consists 
of employing a thrcc-constitucnt target, mih the 
relative abundances of all constituents known. We 
argue that the equilibrium plasma model implies a 
thorough mixing of the plasma constituents and 
that, as a result, target and particle constituents 
arc interchangeable in the equations. Denoting 
the target species by (1), (2), and (3), we treat 
(1) and (2) as “knowns,” using a 
measurement to compute T according to equation 
(4), and (3) as the '"unknown,” computing 
from the temperature calculation and 
equation (5). Of course, the actual value of 
is known; this known value can be 
compared to that computed from equation (5) , 
and good agreement is evidence for the validity 
of the model. W^e noted that the evidence is not 
quite so strong as in the first approach, because 
now the particle properties do not directly enter 
into the calculations and hence the assumption 
of mixing, which is central to this work, remains 
an assumption. 

In spite of the above limitation, it is important 
to observe that the state of knowledge of the 
impact ionization effect with respect to basic 
theoretical considerations is sufficiently primitive 
at the present time that any direct evidence for 
the validity of a proposed theoretical model of the 
process w^ould represent a major advance. This 
wras an important consideration in our decision to 
perform the experiment to be d(^ribcd. 

EXPERIMENTAL APPARATUS AND 
TECHNIQUES 

All of the experimental data on the impact 
ionization effect W’cre acquired with the TRW 
electrostatic microparticle accelerator described 
by Friichtcnicht (1962). This device is a two 
million volt Van de Graaff generator in which the 
high voltage terminal has been modified to accept 
microparticle charging and injection equipment. 
It is capable of charging and accelerating a variety 
of micron size range particulate materials to high 
velocities. The final velocity obtained is dependent 
upon the size of the particular particle, the 
material used, and the accelerating potential. 


Particle velocities range dowm to about 1 km/s 
for all materials. The upper velocity limit is, for 
practical purposes, dependent upon the elec- 
tronics used for particle detection. For carbonyl 
iron particle, the maximum detectable velocity 
is about 50 km/s; for low^er density materials, the 
maximum velocity is much higher. Particles may 
be accelerated either continually or in single 
bursts. 

Particles from the accelerator pass through a 
particle charge-velocity detector which consists 
of a cylindrical drift tube of knowm length and 
capacitance and an associated amplifier. The 
rectangular pulse derived from this detector as 
the charged particle traverses its length provides 
the necessary data for computing particle mass 
and velocity. The details have been described 
elsewhere (Shelton ct al., I960; Hansen and 
Roy, 1966). 

The charged particle ‘‘beam” from the accelera- 
tor has a broad distribution due to particle size 
variation. For experimental convenience the beam 
is further conditioned by means of particle velocity 
(i.e., transit time) selection equipment (Roy and 
Becker, 1971) which may be set to limit the 
particles reaching the experimental area to only 
those having a velocity lying within a narrow 
velocity range. All others are electrically de- 
flected. 

A detailed block diagram of the cxpenmcntal 
apparatus including a schematic representation 
of the target and collector assemblies is shmvn in 
figure 1. The tar^t and collector assemblies are 
housed in a cylindrical vacuum enclosure about 
45 cm in diameter by 125 cm long. Particles from 
the accelerator traverse a velocity selection 
system (not shown in the figure) and then a 
charged particle detector prior to entry into the 
main experimental volume. The particles pass 
through an aperture in the ion collector and 
proceed onward to impact on the target assembly. 

The target assembly is mounted at the down- 
stream end of the system (most distant from the 
accelerator) such that simulated micrometeoroids 
produced by the accelerator will impact near the 
center of the target and normal to its surface. 
Since particles from the accelerator arc confined 
to a 6.3 mm diameter circle, a 5 cm diameter 
target was considered ample for these tests. The 
target is mounted inside a shielded housing as 
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Figure I. —Detailed block diagram of the experimental configuration, with the internal target 
and collector assemblies depicted schematically. Signal outputs from the experiment are 
Illustrated in (b) ; the upper trace is the sum of particle detector, target, and collector signals 
with sweep speed set at 20 /xsec/div; lower trace, target and collector signals at 5 ^c/div. 


shown in the figure, with the extraction and 
acceleration grid mounted approximately 1 cm 
from the target surface. Cable connectors \vith 
hermetic seals provide means for introducing 
the target bias voltage (+3000 V) and with- 
drawing the target signal from the vacuum 
chamber. 

The ion collector mechanical assembly is 
mounted at the opposite (upstream) end of the 
vacuum chamber at a distance of approximately 
1 m from the target. The collector plate is about 
30 cm in diameter and is centered on the axis of 
the particle beam. The entire assembly has a 
small hole (about 1.25 cm) at the center W'hich 
allows the particles to pass. The collector plate is 
surrounded by a grounded shield. In order to 
maintain electrical shielding \vhilc permitting 
ions to enter, the side of the shield w+ich faces the 
target is made from screening cloth of reasonably 
high transmittance. The large diameter collector 
plate is required to collect a substantial fraction 
of the ions originating at the target, due to 
expansion of the beam as a result of the initial 
thermal energy of the ions. Again, a hermetically 


sealed cable connector allows the collector signal 
to be extracted from the vacuum enclosure. 

The photographic record of a typical event is 
shown in figure 1(b). This particular oscillogram 
records the data from the impact of an iron 
particle on a PbZr 03 -PbTi 03 target at 36 km/s. 
The sequence of events begins when a particle 
from the accelerator is found by the time interval 
selection equipment to have a transit time within 
the preselected transit time (velocity) boundaries. 
This event causes the bias voltage on a particle 
deflector to be removed, thereby allowing the 
selected particle to enter the experimental region. 
Just prior to arrival of the particle at the particle 
detector depicted schematically in figure 1, a 
proportional delay generator produces an output 
pulse to initiate an oscilloscope sw'cep that will 
present the signal amplitude and transit time 
information produced by the particle as it passes 
through the detector (upper trace on the oscil- 
logram). This information allow's the particle 
parameters (velocity, mass, and radius) to be 
calculated for each impacting particle. The 
particle then passes through the collector aperture 
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enroute to the target. Just prior to arrival of the 
particle at the target, a second sweep is initiated 
which will present the target impact signal and 
the ion charge and transit time information 
(lower trace of the oscillogram). The signal ob- 
tained from the target is essentially a step func- 
tion for impacts above about 10 km/s, while that 
from the collector rises in staircase fashion. Either 
or both waveforms may be differentiated to 
produce sharp pulses at impact and/or upon 
arrival of each ion group at the collector. 

In Figure 1(b), the lower trace is the sum of 
the target and collector signals at a sweep speed 
of 5 ^sec/division. The initial negative step 
represents charge leaving the target, while the 
following positive-going steps represent the same 
charge arriving at the collector, noiv separated 
into individual charge “packets” distinguishable 
as ions of the constituent elements of particle and 
target material. The upper trace is the sum of the 
signals from the particle detector, the target, and 
the collector. The upper trace is presented at a 
slower sweep speed (20 *isec/division) , insuring 
that the w^aveform being measured is due to the 
impact of the detected particle. 

A total of 33 separate impact events on a lead 
titanate-lcad zirconate target were recorded u.sing 
the apparatus shown in figure 1. Of these, 12 
particles were nickel-aluminidc, seven were 
molybdenum boride, and the remainder w-ere iron. 
The impact velocity range extended from 16.8 to 
46.8 km/s. For each event the total number of 
lead, zirconium, and titanium ions reaching the 
collector was measured. 

DATA ANALYSIS 

All data from impacts of Fe, MoB, and NiAl 
particles onto the PbZrOa/PbTiOs target sarnplc 
w'ere analyzed in accordance wdth the equilibrium 
plasma model discussed in the section on the 
impact model. There we saw that two of the 
target elements were to be treated as knoivns and 
used to compute the relative abundances of the 
third element treated as an unknown. The ioniza- 
tion energies of Pb, Ti, and Zr are, respectively, 
7.42 eV, 6.82 eV, and 6.84 cV. Examination of 
equation (3) shows that the two control species 
(1) and (2) must have an appreciable difference 
AEiU'« in ionization energy in order for this 


equation to yield meaningful results. This clearly 
required that Pb be one of the control elements. 

Ti w'as arbitrarily selected as the second control 
element, meaning that Zr became th(; “unknown.” 
The appropriate partition functions for all 
three elements were available in the literature 
(Aller, 1963) at four fixed temperatures: 3600, 
5700, 8000, and 10 000° K. Using these values and 
the ionization energies of the elements, we were 
able to compute values of the functions 

f[T)=f'{T) cxp( — (6a) 

and 

g(7) =g'{T) exp — Fj‘”)/fcT] (6b) 

where superscripts (1), (2), and (3) now signify 
Pb, Ti, and Zr, at the four given temperatures. 
On the assumption that the partition function 
ratios arc smooth, wnll-bchaved functions, we 
used the computations to plot/(T) and g{T) vs T. 
Equations (3) and (5) can be rewritten 

iV/»/iV.®=-4/(T) (7a) 

and 

= (iV/’’AV/‘>)sf(T) (7b) 

The experiment gave measured values of 
iVi(Pb)A¥<(Ti), which equals 
iv'i(Zr)/iV'i(Pb), which equals 
A is given by A„(Pb)/iSr„(Ti) and is knowm, the 
former value could be used to compute /(T) . From 
the plots, T and then g ( T) were found, after which 
iV«(Zr)/iV„(Pb) could be evaluated. 

The “true” relative abundances of Pb, Ti, and 
Zr in the target sample w’ere determined by 
electron microprobe analyses. To check for 
inhomogeneitics, readings W'erc taken at four 
separate points on the impact surface using an 
incident electron beam diameter on the order of 
one micron, with an analyzed depth also on the 
order of 1 micron. The composition found for the 
four points (table 1) indicates good homogeneity. 
The mean composition of the crystal W’as found to 
be as follow’s: 

Pb = 58.62± .55 atomic percent 
Ti = 21.17zb.l6 atomic percent 
Zr= 19.73 ±.67 atomic percent 

where all quoted error values are 2cr limits. In 
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terms of relative abundances, 

iVn(Ti)/iVn(Pb) =0.361 it .006 (8a) 

and 

Nn{Zr)/Nn(Ph) =0.337±.015 (8b) 

The mechanics of data reduction to A^n(Zr)/ 
A^n(Pb) were straightforward and should be 
obvious in the light of the preceding discussion. 
The graphical method by which /(T) values were 
used to obtain g{T) naturally contains a number 
of approximations whose effects on the accuraFcy 
of the experiment are difficult to determine. 
Existing theory is sufficiently complete to permit 
the use of a high-speed digital computer to 
separately evaluate the partition functions Ui and 
Un for almost any material at any temperature. 
Ideally, such a computation could have been 
performed for a large number of temperatures, 
after which the computer could have evaluated 


f{T) and g(T) at these temperatures. Interpola- 
tions necessary to obtain the functions at any 
temperature could also have been included in the 
program. Such techniques will actually be em- 
ployed to reduce cosmic dust analyzer in the space 
experiment, but they did not appear justified in 
the pr^ent context due to time and cost factors. 
After examining all aspects of the current experi- 
ment, we conclude that our graphical approxi- 


Table l— Electron Microprobe Analysie of Target Sample 
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Figure 2.— Computed impact plasma temperature vs particle impact velocity: various particles, 
PbZrOa-PbTiOs target. Data points shown are the computed result using the known relative 
abundance of two target elements together with the impact ionization model. Computed 
average and median impact temperatures as indicated. 
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mation has not materially affected its significance 
or validity. 

Figure 2 is a plot of the impact plasma tem- 
perature computed from equations (7a) and (4) 
as a function of impact velocity. The data scatter 
widely and there is no observable trend ^^ith 
increasing velocity j NlAl impacts appear to 
produce lower plasma temperatures on the average 
than do Fe impacts. (The model would predict 
such a trend due to the lower ionization potential 
of aluminum.) These results are similar to what 
has been observed with all other particle/ target 
combinations employed during this program, and 
they tend to confirm that mechanical details of 
the materials at the impact sites (yield strengths, 
angle of impact, hardness, etc.), all of which arc 
essentially random variables in the size range of 
the particles impacted, are more important in 
determining the amoimt of energy available for 
material vaporization and ionization than is the 
total energy of impact. 

Remember, however, that the equilibrium 
plasma model says nothing at all regarding the 
relationship, if any, between impact velocity and 
plasma temperature. In fact, when the seemingly 
scattered data are analyzed further to compute 
An(Zr)/iVn(Pb) according to equation (7b), 
surprisingly good agreement is found from one 
point to the next. Table 2 summarizes the velocity, 
plasma temperature, and relative Zr abundance 
iVn(Zr)/iVn(Pb) for all of the data; the average 
Zr abundance is 

<iV„(Zr)/iV„(Pb) > =0.38±.12 (9) 

where the error value is a 2<r limit. Comparison 
with equation (8b) shows very good agreement, 
and the value of the standard deviation (<t = 0.06) 
for the 33 points is reasonably small. In addition, 
no one data point deviates more than about 40 
percent from the average, vrhich is rather remark- 
able given the fact that the data were acquired 
over a considerable period of time with three 
different particle materials and the method of 
analysis was rather crude in several aspects, as 
wc have noted. 

The sensitivity of this system for predicting 
elemental relative abundances from ion ratio 
measurements is illustrated in figure 3, where we 
have plotted the Zr/Pb ion ratio as a function of 
the Ti/Pb ion ratio. The three curves have been 


computed for different assumed compositions of 
the target material, using the four fixed tempera- 
tures at which the partition functions were avail- 
able; the middle curve corresponds to the actual 
composition given by the microprobc analysis, the 
upper curve to a 50 percent higher Zr content, 
and the lower curve to a 50 percent lower Zr 
content. (In both cases, the Ti content has been 
altered in inverse proportion.) The measured ion 
ratios are seen to cluster quite closely to the middle 
curve, and there is little difficulty in determining 
that the Zr content of the substance is much closer 
to 20 percent than to either 10 or 30 percent. 

Figure 4 is a plot of the computed Zr/Pb ratio 
as a function of impact velocity. It is included as 
one possible check on the validity of the model 
and the method of analysis, since if both are 
correct the ratio iVn(Zr)/iVn(Pb) should be 
independent of v. The plot shows this to be the 
case. 

CONCLUSIONS 

The results presented in the preceding section 
provide fairly convincing proof that the use of an 
assumed thermal equilibrium plasma model can 
be used to process data from an impact ionization 
time-of-flight mass spectrometer in order to 
convert the raw' ion data to relative abundances 
of the elemental constituents of cosmic dust 
particles. We do not wish to imply that the 
equilibrium condition is actually met; even in- 
tuitive arguments w^ould suggest that a plasma 
cannot be in thermal equilibrium shortly before 
the plasma becomes collisionless, because the 
collision frequency must approach zero in some 
asymptotic fashion. Hovrever, a * ‘quasi-equilib- 
rium” plasma is conceivable, and the anatysis 
wnuld then yield an “effective” plasma tempera- 
ture. On the basis of the data presented in figure 
3, the model appears to be relatively sensitive in 
computing relative abundances. Considering the 
application and the precision requirements that 
can sensibly be placed on this kind of flight 
experiment, the rudimentary model appears to be 
quite satisfactory. 

As mentioned earlier, a remaining question is 
verification that the particle material becomes 
thoroughly mixed with the target material, which 
is a requirement if a common plasma temperature 
is to be assigned to both materials. Our search for 
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Table 2. — Summary of Experimental Data 


Date and 
particle material 

Temperature 

CK) 

Velocity 

(km/s) 

iV„(Zr)/JV„(Pb) 

Squared deviation 

9/17/70 

4480 

29.4 

0.255 

0.01440 

MoB particles 

4200 

24.8 

0.354 

0.00044 


11 500 

27.3 

0.374 

0.00000 


8050 

24.0 

0.408 

0.00109 


4050 

33.8 

0.375 

0.00000 


3300 

28.9 

0.250 

0.01562 


4270 

28.4 

0.265 

0,01210 

9/22/70 

5550 

29.8 

0.423 

0.00230 

NiAl particles 

4350 

23.4 

0.295 

0.00640 


6900 

19.4 

0.239 

0.01850 


5150 

27.5 

0.419 

0,00194 


3925 

34.1 

0.371 

0.00002 


4600 

29.3 

0.393 

0.00032 


4800 

31.8 

0.403 

0,00078 


4630 

37.0 

0.374 

0.00000 


4800 

29.6 

0.444 

0.00476 


6050 

25.5 

0.323 

0.00270 


4900 

28.3 

0.384 

0.00008 


3920 

25.9 

0.357 

! 0.00032 

1/16/71 

4850 

19.6 

0.314 

0.00372 

Fe particles 

5100 

18.5 

0.373 

0.00000 


7650 

18.1 

0.419 

0.00194 


5300 

16.8 

0.375 

0.00000 


4800 

17.2 

0.391 

0.00026 

3/4/71 

7650 

39.5 

0.534 

0.02528 

Fe particles 

7280 

26.3 

0.401 

0.00068 


8180 

46.8 

0.394 

0.00036 


7950 

29.8 

0.391 

0.00026 


6470 

29.7 

0.367 

0.00006 


7200 

30.0 

0.432 

0.00325 


8970 

31.4 

0.449 

0.00548 


10 870 

30.7 

0.430 

0.00302 


5750 

35.7 

0.406 

0.00096 




Av=0.375 

Sum =0.12705 


Figure 4— Computed zirconium/lead relative abundance .V n(Zr)AV„(Pb) vs particle impact 
velocity. Zirconium/lead relative abundance is computed using measured experimental values 
and the impact ionization model. Each computed value for relative abundance is plotted 
against the particle impact velocity from which the experimental values were obtained. 
Target surface analysis, as determined by electron microprobe, is as indicated to the right of 
the figure. 
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Figure 3. — Zirconium /lead ion ratio A i(Zr)/A^i(Pb) vs titanium /lead ion ratio A i(Ti)/A^t(Pb). 
Data points are derived from the ratio of the indicated ion groups, the amplitudes of which 
are measured from the collector waveform. Each point represents a single impact event. The 
curves shown are computed using the indicated target composition for each and the equations 
which characterize the impact ionization model. Bending of the curve is caused by the tem- 
perature dependence of the partition functions for ions and neutrals of each of the subject ele- 
ments. 
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satisfactory two-element particle materials will 
continue in an effort to answer this question. 

Another practical task that must be accom- 
plished is the selection of the control elements in 
the target. The control elements can, in principle, 
take several forms; a substrate can by some means 
be impregnated with the two control elements, a 
binary compound material can be used, or vapor 
(or elect roly tically) deposited films can be em- 
ployed. Whatever the form, homogeneity of the 
materials on a microscopic scale is a necessary 
condition. Several other factors must also be 
considered: 

(1) The ionization potentials of the two 
elements must differ by about one electron volt in 
order to obtain a substantial change in the ratio 
of the ion signals from the control elements as the 
plasma temperature changes. 


(2) The ionization potentials must be low 
enough to ensure a sufficient degree of ionization 
to obtain detectable signals. 

(3) The mass of the atoms should lie outside of 
the mass range of elements expected from the 
cosmic dust particles. 

(4) The materials must be stable when exposed 
to the space environment. 

In conclusion we would suggest that laboratory 
experimentation for the purposes of f urther 
verification (and possible improvement) of the 
model and selection and test of target materials be 
continued. An early flight test of even a simplified 
version of a cosmic dust analyzer would be a 
highly desirable undertaking. Such a flight test 
would provide a measure of the signal magnitudes 
from naturally occurring cosmic dust particles 
and could crudely identify the major elemental 
constituents. 
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